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Abstract. The theory of highly super-Eddington black hole accretion, M S> MEdd, 
was developed in the 1980s in Warsaw by Paczyhski and his collaborators in 
terms of Polish doughnuts 1 , i.e. low viscosity, rotating accretion flows that are opti- 
cally thick, radiation pressure supported, cooled by advection, and radiatively very 
inefficient. Polish doughnuts resemble fat tori with two narrow funnels along rotation 
axis. The funnels collimate radiation into beams with highly super-Eddington lumi- 
nosities. Slim disks 1 introduced later by Abramowicz, Lasota and collaborators first 
m [TU). and more fully in have only moderately super-Eddington accretion rates, 
M > AfEdd, rather disk-like shapes, and almost thermal spectra |12 | -|14 | . 

1 Introduction 

Is there an upper limit to the accretion rate M, preventing astrophysical black 
holes to swallow all the matter infalling into them? Many astrophysicists are 
convinced that the answer to this fundamental question should be: yes, the limit 
exist and it is given by the standard Eddington accretion rate, 

M Edd = ^ = 1.5x 10 17 (J£) [g sec" 1 ]. (1) 

However, because there is a clear observational evidence for super-Eddington lu- 
minosities, e.g. from X-ray binaries (SS433, GRS 1915+105), and because super- 
Eddington luminosities powered by accretion imply highly super-Eddington ac- 
cretion rates, this almost unanimous consensus 2 is to me rather befuddling. 

I review here fundamental theoretical predictions about super-Eddington ac- 
cretion: (1) super-Eddington luminosities are typical for rotating, radiation pres- 
sure supported black hole accretion flows, (2) these flows have very small effi- 
ciency and therefore they must have highly super-Eddington accretion rates, 
(3) super-Eddington accretion does not necessarily imply strong outflows. 

These predictions are solidly based on standard physics, and it should not be 
surprising that recent 3D hydro and MHD supercomputer simulations provide 
many beautiful illustrations for them |15 | -|18 | . 

1 The name Polish doughnuts was coined by M.Rees, and slim disks by A.Kolakowska. 

2 All works on cosmic structure formation reported at this Conference, used MEdd as 
the upper limit for the black hole growth rate. 
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1.1 Sub-Eddington accretion: standard thin disks and adafs 

The sub-Eddington accretion is far more familiar to astrophysicists. Well known 
theory predicts for M <§C MEdd two types of stable accretion, depending on the 
optical depth, r. 

When r ^> 1, accretion is described in terms of the famous standard thin 
disk model by Shakura & Sunyaev JH]- Standard thin disks are supported by 
gas pressure, cooled by radiation, and very efficient. They are geometrically thin 
in the vertical direction, H/R <§; 1. 

When t <C 1, accretion is described by the adaf 3 model, propheted in 1977 
by Ichimaru in a paper |2()j that was ignored by everybody for twenty years 4 . 
Only after their rediscovery in the mid 1990 by Narayan 5 & Yi [22], and 
by Abramowicz, Chen, Kato, Lasota & Regev [21], adafs started to be intensely 
studied by many authors (see reviews in [25], [2H])- Adafs are cooled by advection. 
They are very inefficient, geometrically thick, H/R ~ 1, and very hot (close to 
the virial temperature). Because of their very low efficiency, they are much less 
luminous than the standard thin disks. 

2 Fundamental limits for power: Planck and Eddington 

The Planck power (i.e. power expressed in Planck's units) equals 

ipianck = ^ « 10 58 [ergsec- 1 ] = 10 52 [Watts]. (2) 

Rather surprisingly, it does not depend on the Planck constant h. The maximal 
energy available from an object with the mass M (and gravitational radius Rq = 
GM I c 2 ) is -Emax = Mc 2 . The minimal time in which this energy may be liberated 
is imin = Rg/c- Thus, the maximal power L max = E max /t min = c 5 /G = L Planck . 
This is the absolute upper limit for power of anything in the Universe: all objects, 
phenomena, explosions, and evil empires 6 . 

For example, a sphere with radius R containing blackbody radiation at tem- 
perature T radiates power L <~ (caT A )R 2 . The gravitational mass of the radiation 
inside the sphere is M ~ (aT 4 /c 2 )R 3 , where a = 87r 5 fc 4 /15c 3 ft, 3 is the radiation 
constant. Thus R G - (aT 4 G/c 4 )R 3 , and from R G < R one gets L < c 5 /G. 

Consider a stationary object with mass M in which gravity and radiation are 
in equilibrium. Sikora 27, noticed that the upper limit for the object radiative 
power may be expressed by LEdd = ^PianckJC'grav/.E'rad- Here JCgrav is the object 

3 The name adaf was coined by J.-P.Lasota. 

4 However, many elements of the adafs theory were present in the influential "ion tori" 
paper by Rees, Phinney, Begelman & Blandford |21) . 

5 R. Narayan did more than anybody else to develop the adaf idea into a useful and 
mature astrophysical theory. 

6 A power needed for the Creation was the rate at which the Big Bang transferred 
energy from a pre-Planck to the post-Planck state. For the reason outlined here, no 
more than 10 52 Watts was needed to create the Universe. 
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total effective gravity cross section, and 27 ra d is its total radiative cross section. 
Sikora's expression is the most general version of the Eddington limit. Eddington 
himself, first in [2H], and then in considered a much more specific case 

of a radiation pressure supported star 7 , assuming that radiation interacts with 
matter by electron scattering (thus S Ta d = NaT in Sikora's expression) , and that 
effective gravity is provided by the Newton gravity alone (thus S STa ,v = 4ttRq). 
Here N is the number of electrons in the object (TV = M/mp for pure hydrogen, 
mp is the proton mass), o~t = 87re 4 /3rn^c 2 is the Thomson cross section (e is the 
electron charge and m e its mass). In this specific case, Sikora's elegant argument 
immediately gives the standard formula for the Eddington limit, 



i E dd = ipianck^ 1 = — = 1.4 x 10 38 — [erg sec" 1 ] . (3) 

^rad CTT \Mq/ 

The Eddington limit for a spherical, non-rotating, homogenous object was dis- 
cussed from a modern perspective by Joss, Salpeter & Ostriker 021- I n addition 
to the discussion given there, let us note that, obviously, the limit for radiative 
power may increase above the standard Eddington limit, in objects with a grav- 
itational cross section greater than the standard one, or with a radiative cross 
section smaller than the standard one 8 , 

<^grav > ^grav = ^Rq, Kud < ^rad = — CTT ■ (4) 

77lp 



In 1913, shortly before the war, Bialobrzeski 33 pointed out that radiation pressure 
may be important in stellar equilibria (see also |34 ) -|3(i ) 'l. Eddington independently 
discovered the same three years later |28l . already during the war. Eddington, as a 
Quaker, avoided an active war service and was able to continue his research during 
the war years 1914-1918 at Cambridge, as the Director of the Observatory there. Only 
after the war ended, the two scientists could communicate, and Eddington wrote to 
Bialobrzeski: / congratulate you on having been apparently the first to point out the 
large share of radiation pressure in internal equilibrium of stars (quotation after |37|h 
Bialobrzeski's contribution (and priority) was remembered by astronomers still in the 
1930s, (see the 1936 monograph by Tiercy |38| and the 1939 book by Chandrasekhar 
|39|h but today is almost totally forgotten — Chandrasekhar |4U| makes no mention 
of it in his 1983 book on Eddington. ( A.K. Wroblewski ,41, provided most of historical 
information for this footnote.) 

In no circumstances LBdd can grow above Lpianck- I noticed I I 31 that the two limits 
are equal when M ta Nomp, with No = ^Vn irac /3, where the Dirac number Afcirac = 
e 2 /Gm e mp equals the ratio of Coulomb's to Newton's force between electron and 
proton. Then iV » iV E dd = 136 x 2 256 w 1.6 x 10 79 , where JV Ed d is the Eddington 
number that played an important role in Eddington's Fundamental Theory |44| . The 
number was introduced by Eddington's immortal statement, / believe there are 15 
747 724 1 36 275 002 577 605 653 961 181 555 468 044 717 914 527 116 709 366 231 
425 076 185 631 031 296 protons in the universe and the same number of electrons. 
One thus may write, L-^dd = (N/N-Edd) ipianck- A very Eddingtonish connection 
indeed, embracing his luminosity and his number, of which he was not aware. 



4 M.A. Abramowicz 



Several types of astrophysical sources give a clear observational evidence for 
super-Eddington luminosities. I will mention here only two specific types that 
provide examples for the two possibilities listed in (@J. 

Novae outbursts. They are due not to accretion but to thermonuclear power 
|45| . Luminosity increases within a few hours by factors of 10 4 , and stay for a very 
long duration at a clearly super-Eddington level, L > lOLEdd- Shaviv 05] > SHI 
shows that the observed increase of the radiative power over the Eddington limit, 
may be attributed to I7* ad < S la d, because in a locally inhomogeneous medium, 
the ratio of averaged radiation force to emitted flux, F ra d/ / ra d = c* ad /c < o~t/c. 

X-ray binaries. Some of them show super-Eddington luminosities, powered, 
most likely, by super-Eddington accretion: in high-mass binaries undergoing a 
thermal-time-scale mass transfer, e.g. in SS433, |47| . and in the low-mass binaries 
during long-lasting transient outbursts, e.g. in GRS 1915+105, |4"5], 

The rest of my review is devoted to the discussion of a possibility that the 
reason for super-Eddington luminosities in the X-ray binaries and similar objects 
may be the fast differential rotation of accretion disks present in these sources. 
It strengthens the effective gravity, so that S* rav > AttRq. 

3 Eddington limit for rotating bodies 

I now consider the Eddington limit for a perfect-fluid rotating body in equilib- 
rium, following the line of arguments first presented by Abramowicz, Calvani & 
Nobili jSJ. Having in mind most general astrophysical applications, I will consider 
two topologically different cases: 

Rotating stars. The surface of the star has a topology of the sphere. The 
whole mass M is included in the sphere. 

Accretion disks. The surface of the disk has a topology of the torus. The mass 
of the disk Md; s k is contained in the torus, but the mass M ccntr of the central 
black hole is outside. The total mass M — M ccntT + Mdisk ~ M CC ntr, because 
-Mccntr 3> -Mdisk- In accretion disk theory it is customary to neglect the mass of 
the disk, so formally M = M ce ntr, and Mdisk = 0. The Eddington limit always 
refers to the total mass M. 

Let f ra d be the local flux of radiation somewhere at the surface of the body, S. 
The corresponding radiative force is F ra a = (c ra d/c)frad- Let F cr f = F grav + F rot , 
be the effective gravity force, with F grav = mV# being the gravitational force (<P 
is the gravitational potential), and with F rot = m(f2 2 r) e r being the centrifugal 
force (J? is the angular velocity, r is the distance from the axis of rotation, and 
e r = Vr a unit vector showing the off-axis direction). 

The necessary condition for equilibrium is F ra d < F e ff. From this one deduces 
the Eddington limit for rotating perfect-fluid bodies, 

L = f f rad • dS = — f F rad ■ dS < — f Fes ■ dS = Lgk (5) 

JS a rad Js °rad Js 

Using the Gauss theorem to transform the surface integral of F c ff into a volume 
integral of V • F c g , Poisson equation V 2 <£ = AnGp , and introducing the specific 
angular momentum i — fir , one gets after twenty or so lines of simple algebra, 
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^Edd — ^Edd [ ^Lss + -^shear _ ^vorticity ] ) (6) 

where the dimensionless number X^ ass depends on whether the body is a star, 
or an accretion disk, 

^ j. f 1 for stars, 

*Lss = m / pdV = I (7) 
Jy [ for accretion disks, 

and where -^ 2 hear , Vorticity are dimensionless, necessarily positive quantities, 
representing shear and vorticity integrated over the whole volume of the body, 

^=l^//(™-™W (8) 

X ^ = ^j v r-\m.m )d v. ( 9) 

Shear increases the Eddington limit, and vorticity decreases it. 

The rotation of astrophysical objects is far from simple, but an insight could 
be gained by considering a simple power law for the angular momentum dis- 
tribution, £(r,z) — £§r a '■> with £q and a constant. Rigid rotation has a = 2, 
Keplerian rotation a = 1/2, and constant angular momentum rotation a = 0. 
K is ^shcarAvorticity = (<* ~ 3) 2 / 'a 2 . This means that A 2 hcar > A 2 orticity when 
a < 3/2. 

Rotating stars have A 2 hcar < A 2 orticity , because they rotate almost rigidly. 
Thus, for rotating, radiation pressure supported stars, L w i^dd < ^Edd always. 
Contrary to this, constant angular momentum tori arc dominated by shear, 
^shcar ^ 1 ^ ^vorticity an d consequently, when they are radiation pressure 
supported, L w Lgjj d > L E dd- 



4 Polish doughnuts 

For a simplicity of presentation, let me assume that £ — £o — const, i.e. that the 
angular momentum is constant in the whole body. One does not know a priori 
what is the actual distribution of the angular momentum inside the body, as 
this depends on the nature of viscosity. Paczyhski explained why a physically 
realistic distribution of angular momentum must indeed be close to £ = const 
in the inner part of the flow, and why most likely it approaches the Keplerian 
distribution, £{r, z) = £k(t) = (GMr) 1 / 2 far away from the center. 

Adopting Paczyhski's argument, one may assume that the local physical 
properties of the innermost part of black hole accretion flows are rather well 
described by the model with £ = const, but one must be careful with the physical 
interpretation of the £ = const assumption at large radii. I will return to this 
point later in this section. 

Assuming constant angular momentum, and using Newton's expression for 
gravitational potential <S> = — GM/(r 2 + z 2 ) 1//2 in cylindrical coordinates [r,z,(p], 



6 M.A. Abramowicz 



I deduce the shape of equipressure surfaces P = const from the Bernoulli equa- 
tion, 

CM 1 f 2 

uivj_ i < ^ (p) = congt (1Q) 

(r 2 + z 2)2 2r 2 

Here W{P) — — J dP/p. This equation cannot be obeyed at the rotation axis 
(where z / 0, r = 0), which has an obvious, but important, consequence: no 
equipressure surface cross the axis. For a constant angular momentum fluid, 
equipressure surfaces must be toroidal, or open. The marginally open surface 
has just one point (r = oo, z = 0) at infinity. This particular surface encloses the 
largest possible torus. From (fTUfl it is obvious, that in this case W(P) — W(0) — 
0. 

Maximal pressure locates at a circle z — 0, r = r$ = GM/£q. Using the radius 
r as a scale, £ = r/r , r\ = z/r , w — W/(GM/r ), one may write equation ltTU|l 
in the dimensionless form, and solve for rj = 7j(£; w) to obtain the explicit form 
of all equipressure surfaces, 

?? 2 = Q 2 (0=4£ 4 (l-2£ 2 W )- 2 -£ 2 , - 1/2 < w < w s . (11) 

The value w = —1/2 gives the location of the center, and w = ws < the 
location of the surface. For the slender torus ws ~ —1/2, and the fat torus 
ws ~ 0. 




Fig. 1. Contours of equipotential surfaces on the meridional cross section of a constant 
angular momentum torus. Left: a slender torus, ws ~ —1/2. The contours approach 
circles around the locus of the maximum pressure, at £ = 1. Right: a fat torus w s ~ 0. 
The contours tend to concentric circles with the center at £ = 0. Close to axis they 
change into a pair of conical funnels. 
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4.1 The slender torus 

One may introduce toroidal coordinates centered at the circle of the maximal 
pressure by a coordinate transformation, x — [rj 2 +^ 2 ] 1 ^ 2 , a = tan -1 (ji /£,*). Here 
= 1 — £. It is then not difficult to show that in the slender torus limit, i.e. 
when to + 1/2 = e <C 1, it is w = w(x) + 2 (e). This means that at the meridional 
section of the torus, equipressure surfaces are concentric circles, as first noticed 
by Paczyhski 49 . This additional symmetry was used by Blaes jSU] to obtain 
a general analytic solution for all possible slender torus oscillation modes. In 
particular, Blaes demonstrated that there exists a set of non-axisymmetric (i.e. 
with to 7^ 0) global modes SX mn = F mn {x, a) ex.p(—iLu mn t + m<p) of the slender 
torus oscillations, with the eigenfrequencies, 



Because 3(w mn ) < these oscillations are unstable. The growth rate of the 
instability is ~ rnQx- Here Ok is the Keplerian frequency at the torus center, 
and thus the instability is a dynamical one. Indeed, this is the famous instability, 
discovered in the seminal papers by Papaloizou and Pringle |51| . |f>2| . 

4.2 The fat torus 

It is easy to understand shapes of equipotential surfaces in a fat torus (shown in 
Figure Very far from the axis of rotation one has \2w£\ ^> 1. Inserting this 
into one gets rj 2 +^ 2 = 1/w 2 which means that far from the rotation axis the 
equipotential surfaces are spheres with radius and that the outer radius 

of the torus is at r out = ro/\ws\- Spherical cquipotentials are in accord with the 
fact that very far from the axis, £ 3> 1, the centrifugal force ~ ^o/£ 3 is negligible 
with respect to the gravitational force ~ GM/£ 2 . Therefore, the effective gravity 
is determined by Newton's attraction alone, as for spherical stars. For the same 
reason, the radiation power from this part of the surface of a fat, radiation 
pressure supported torus (i.e. Polish doughnut) is one Eddington luminosity, the 
same as from a spherical non-rotating, radiation pressure supported star 9 . 

Note, however, that the asymptotically spherical shape of a fat torus is a 
direct consequence of the assumption £(r, z) = const, which was made ad hoc. If 

9 M. Rees told me in 1980 that this implied that Polish doughnuts, with at least 
one Eddington luminosity from whatever direction, were obviously too luminous to 
describe very "dim" active galactic nuclei, such as radio galaxies: "while apparently 
supplying tremendous power to their extended radio-emitting regions, the nuclei of 
most radio galaxies emit little detectable radiation." He and his collaborators at 
Cambridge found later a possible solution to this puzzle in terms of the ion pressure 
supported tori |21| . The ion tori have the same shapes as the Polish doughnuts, in 
particular funnels, but much lower, indeed very sub-Eddington, luminosities. The 
power in jets comes from tapping rotational energy of the central black hole by the 
Blandford-Znajek mechanism |54|. 




(12) 
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one adopts a more physically realistic assumption that asymptotically £(r, z) = 
Iki one may use the standard Shakura-Sunyaev model in its radiation pressure 
version, to get the asymptotic shape of the fat torus, 

Zoo = JZ^M. (13) 



Closer to the axis, \2w£\ ~ 1, which means that £ 2 ~ —l/2w and this together 
with gives rj 2 /^ 2 = (l/\w\) - 1 = (l/sin 2 0) - 1, i.e. that closer to the 
axis, equipotential surfaces corresponding to w ~ ws ~ have conical shapes 
with the half opening 9 ~ y/w. The surfaces are highly non spherical because 
centrifugal force dominates. Integrating effective gravity along the conical funnel 
is elementary, and one gets [3] that L^ d /LEdd = (1/2) ln(l/|ws|). This estimate 
may be used to find the total luminosity for the radiation pressure supported 
fat torus, i.e. a Polish doughnut, 



In 1^ = i.i5 i og 1^1 . (14) 
^Edd ^Edd 2 \ r J \ r J 

It should be clear from our derivation, that the logarithmic scaling of the lu- 
minosity with the torus size is a genuine property of the Polish doughnuts, 
including those that have a non-constant angular momentum distribution. The 
logarithm in l|14[) is of a crucial importance, as it prevents astrophysically re- 
alistic doughnuts (i.e. with r out /ro < 10 6 , say) to have highly super-Eddington 
luminosities. Thus, the theory predicts that for such "realistic" fat tori, only a 
slightly super-Eddington total (isotropic) luminosities, Ltotal > 7LEdd, may be 
expected. 

However, because the funnels have solid angles 9 2 ~ ro/r ut, radiation in 
the funnels may be, in principle, collimated to highly super-Eddington values 
L co \\l '^Edd = ~ foutAo ^ 1- This simple estimate agrees with a more detailed 
modelling of the Polish doughnuts radiation field by Sikora [55] and Madau [SSj 
who obtained < 10 2 for disks with r out /ro ~ 10 2 . A typical value that follows 
from observational estimates for non-blazar active galactic nuclei, e.g. by Czerny 
and Elvis |S7] and Malkan [SHj, is ~ 10, but of course for blazars and other 
similar sources (e.g. for ULXs, if they are powered by stellar mass black holes, 
as argued by King [SHj), it must be » 10. 

Such high values of are also consistent with the idea, suggested by Paczyh- 
ski [HO] and independently by Lynden-Bell [HI] , that relativistic electron-positron 
e~e + jets may be very effectively accelerated by the radiation pressure in the 
fat tori funnels 10 . Note, that if the flux in the funnel is times the Eddington 
flux, the e~e + plasma feels the "effective" radiative force corresponding to the 
Eddington ratio m p /m e « 10 3 times greater 11 . Detailed calculations (e.g. \<52\ . 
|55| . demonstrated that indeed the e~e + jets may be accelerated in funnels 

10 Lynden-Bell called this an "entropy fountain" . 

11 Note that in Shaviv's explanation of the situation during nova outburst, the luminos- 
ity is physically very close to the Eddington one (in the sense of Sikora's definition), 
although its value is greater than the standard Eddington limit. For the e~e + jets 
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up to the Lorentz factor 7 < 5. However, if jets are initially pre-accelerated 
by some black-hole clcctrodynamical processes (such as the Blandford-Znajek 
mechanism [HI]) to highly relativistic velocities 7 > 10 6 , they will be decelerated 
in the funnels by the Compton drag, reaching the asymptotic Lorentz factor jri4| . 



Observations tell that 700 < 10 2 , and thus equation (|T5|) suggests that O < 10 3 . 
4.3 Rise and Fall of the Polish Doughnuts 

In the time of their debut, Polish doughnuts could theoretically confirm the ob- 
served super-Eddington luminosities, highly collimated beams of radiation, and 
perhaps even the relativistic speeds of jets, fulfilling another principle attributed 
to Eddington |65| : one should never believe any experiment until it has been 
confirmed by theory. These virtues attracted initially some interest of the astro- 
physical community, but the interest has quickly drained with the discovery of 
the Papaloizou-Pringle instability. It was thought that Polish doughnuts must 
necessarily suffer from the instability and thus, in reality, they cannot exist. The 
important discovery by Blaes [HH] that the Roche lobe overflow stabilizes Polish 
doughnuts against the Papaloizou-Pringle instability, came too late — in the 
advent of numerical simulations of black hole accretion flows. Too late, because 
numerical simulations rediscovered and absorbed many of the Polish doughnuts 
results. Today, these results exist in the consciousness of many astrophysicists 
as a set of several numerically established, important but unrelated facts. They 
do not form a consistent scheme that the Polish doughnuts once offered: clear, 
simple, following directly from the black hole physics 12 . 

In the rest of my review, I recall the most fundamental of these results, almost 
totally forgotten today: the relativistic Roche lobe overflow mechanism. The 
mechanism not only stabilizes Polish doughnuts against the Papaloizou-Pringle 
and other instabilities. It also assures that highly super-Eddington accretion 
rates always imply a very low efficiency of accretion 13 . 

the situation is very much different (opposite): the luminosity is physically very 
super-Eddington, although it may be just slightly above the standard limit. 

12 There are at least three brilliant and very pedagogical expositions of the Polish 
doughnuts scheme: two by Paczyhski himself, |(i7) . |fi8) . and one in the text book by 
Frank, King and Raine |69|. 

13 Soltan's well known argument |7U| shows that on average the AGN efficiency cannot 
be small. The argument does not exclude a possibility of a flip-flop sequence of 
periods of low and high efficiency, rather like in |71| for thin disks. Although the 
possibility of a flip-flop behavior with timescale ~ 10 4 yrs (for AGN) has some 
observational back-up 72 , it was never studied sufficiently deeply. 




(15) 
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5 Efficiency of black hole accretion 

I like to illustrate the physical meaning of accretion disk efficiency in terms of the 
Bekenstein engine 14 . Bekenstein [JI; discussed a black hole engine that converts 
mass to energy with a perfect efficiency, following an earlier unpublished remark 
by Geroch |75j . The engine works by slowly lowering on a strong wire a mass 
m into a black hole with the Schwarzschild radius Rs — 2 Ra — 2GM/c 2 and 
mass M. 

As the mass is lowered to a radius R, the energy E{R) measured at infinity 
is gravitationally redshifted and thus goes down relative to the initial energy 
Eq = mc 2 . The change in energy equals to the mechanical work done by the 
wire back at the engine. The efficiency e = [E — E(R)]/E may be calculated 
from the gravitational redshift formula, e = 1 — (1 — Rs/R) 1 ^ 2 - If the mass could 
be lowered to the horizon R — > Rs, the efficiency would go to 1. Of course, 
the mass cannot touch the horizon, because the tension in the wire would be 
then infinite, and even the strongest wire would break. Real wires, that may 
sustain only a finite tension (see e.g. Gibbons jTCj), would break a finite distance 
R = Ri„ from the horizon. In a region below the breaking radius, R < i?i n , 



the mass falls down freely, keeping its energy measured at infinity unchanged. 
Therefore, the efficiency is determined by the radius R = R ln at which the wire 




Fig. 2. Left Figure shows that angular momentum in the fat torus crosses twice the Ke- 
plerian angular momentum. Right Figure, shows schematically that if the the doughnut 
slightly overflows the Roche lobe, a dynamical mass loss will occur. 

With centrifugal support playing a role of the rope, accretion's efficiency is 
also determined by the radius R = i?j n at which the centrifugal support breaks. 

14 Here I borrow a short description of the engine given in |73) . 
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Consider a collection of free particles on Keplerian circular orbits around a 
black hole. Small viscosity will slowly remove angular momentum, so particles 
will slowly drift inward. However, in the region R < i? ms , there are no stable 
circular orbits, and therefore the centrifugal support breaks: the particles free- 
fall into the black hole. Thus, similarly as in the case of the Bekenstein engine, 
the Keplerian binding energy of particles at the location of the marginally stable 
orbit, R ms (also called ISCO), determines the efficiency of accretion. 

5.1 The relativistic Roche lobe overflow 

The location of R ms is determined by the minimum of the Keplerian angular 
momentum £ = Ik, shown in Figure [SJ left. For a constant angular momentum 
fluid torus, £ — £q, the combined pressure and centrifugal support breaks at the 
radius R m < R ms , defined by £k = £o and shown in the same Figure, by the dot 
left of i? ms on the angular momentum distribution. 

At this location, as discovered by Abramowicz, Jaroszyhski and Sikora 1 , 
one of the equipotential surfaces, called the Roche lobe, self-crosses as shown 
in Figure |21 (right) . It should be obvious that if the fluid distribution overflows 
the Roche lobe, a dynamical mass loss must occur: at the the circle R = i?i n 
the centrifugal and pressure support breaks down, and fluid if free falling in the 
region R < i?i n . 

Again, the efficiency is determined by the Keplerian binding energy at break- 
ing point, i?i n , called the inner edge of the accretion disk. This name is the source 
of a confusion, as some wrongly imagine that the name implies that at the in- 
ner "edge" velocity, density, and pressure must experience a jump, or a sudden 
change. Of course they do not. What does change in the region close to R[ n is the 
nature of the flow: from sub-sonic to super-sonic. Note also that for low angular 
momentum flows, with £ < £k everywhere, there is no Roche lobe, and there- 
fore the inner edge cannot be sensibly defined. Such flows differ considerably 
from the Polish doughnuts, see e.g. [77]. Typically, they assume high viscosity, 
while Polish doughnuts assume low viscosity, less than about 0.03 in terms of 
the Shakura-Sunyaev a-parameter. 

It was shown that the inner edge must locate between the marginally 
stable, and marginally bound orbits, 

fimb < R\n < Rms- (16) 

The marginally bound orbit at i? m b has the same binding energy at infinity, 
W(oo) — W(R m h). This has two important consequences: the efficiency of a 
torus with i?i n = R m h is zero, and such a torus must have its outer radius 
at a very large distance, formally at infinity, i? ut = oo. We have previously 
shown that the total luminosity of a fat torus increases logarithmically with the 
outer radius. This, combined with the accretion rate being proportional to the 
size (fT3")l. yields the logarithmic behavior Ltotal/^Edd ~ hi(M/MEdd)- Indeed, 
Paczyhski 0] found that all fat disks calculated by him obey 
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Fig. 3. Left Figure shows a snap-shot of the velocity field inside a time-dependent accre- 
tion black hole flow. The nature of the flow clearly changes at the inner edge (although 
physical characteristics do not experience a jump there). Far away from the inner edge, 
the flow pattern is very complicated, showing back flows, vortexes, and circulations in 
convective cells. Close to the inner edge, the Roche lobe overflow mechanism organizes 
the flow into a highly ordered, almost radial pattern. Right Figure shows that despite 
all complications, the Roche lobe overflow accretion rate self-regulates according to the 
simple analytic formula 1181 . The points are calculated in supercomputer simulations 
|15| of 3D, viscous, time dependent, black hole accretion flows with different parameters 
and different boundary conditions. The straight line is given by the formula 118H . One 
should appreciate the nearly perfect agreement between this simple analytic formula 
and the most sophisticated, indeed today's state-of-art, numerical simulations. Both 
figures are taken from |15|. 



itotal „ al /l6M 

= 4.8 log -. 

^Edd y M Edd 

This shows that the super-critical accretion is qualitatively different from the 
standard, sub-Eddington, thin disk accretion. For the standard thin disks, the 
total luminosity is directly proportional to the accretion rate Ltotal/^Edd = 
e(r ms )M/M Edd , e(r ms ) = const. 

The Roche lobe overflow mechanism self-regulates the accretion rate, which 
for a polytropic equation of state, P — K p 1+1 / n equals, 

M = A{n)K- n — %g A AW\ n+ \ AW = W s -W in (18) 

with an analytic expression for A(n) explicitly known in terms of the Euler 
gamma function, r(n). Because n = 3 for a radiation pressure supported gas, 
the self- regulation imposed by (jT5|) is a very strong one, M ~ |Z\VF| 4 . Suppose, 
that in the region close to the inner edge, there is a fluctuation in the thermal 
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balance, causing overheating. This will induce an expansion, and an increase in 
the Roche-lobe overflow, which in turn will increase the accretion rate, i.e. the 
mass loss from the region. The heat contained in the mass that is lost will cool 
down the region, thus assuring thermal stability. I found that this mechanism 
of advective cooling caused by the Roche lobe overflow, always stabilizes the 
innermost region of any accretion disk (thin, slim, thick, adaf) against thermal 
and viscous instability. Roughly speaking, the local instabilities have no chance 
to grow in the innermost region of accretion disks, because they are quickly 
washed away by advection caused by the Roche lobe overflow. 
Stabilization of the global Papaloizou-Pringle instability by the Roche lobe over- 
flow, found by Blaes |SU], has a different physical reason: almost a perfect re- 
flection of the Papaloizou-Pringle mode at the inner disk edge is a necessary 
ingredient of the instability. The reflection is not possible in the fast (transonic) 
flow induced by the Roche lobe overflow. Thus, the non-realistic, non-accreting 
Polish doughnuts indeed always suffer from the Paploizou-Pingle instability. In a 
more realistic model of a Polish doughnut, that includes the Roche lobe overflow, 
the Paploizou-Pringle instability does not operate. 



6 Discussion 



6.1 Paczyriski's pseudo-potential 

An accurate and elegant model for a non-rotating black hole's gravity was intro- 
duced by Paczyhski jS] it terms of a Newtonian gravitational "pseudo-potential" , 

, , GM 2GM , . 

4>pw(R) = -— — — , Rs = —jr. (19) 
K — its c z 

The Paczyhski model is remarkably successful: numerous authors used it in their 
calculations of black hole accretion flows. This clever idea cannot be, however, 
applied for: (a) rotating black holes, because of the Lense-Thirring effect 15 , and 
(b) self-gravitating fluids, as X/ 2 <Ppw ^ 0. 



6.2 Standard thin disks are inconsistent with M > Msdd 

Properties of accretion flows with M > MEdd are sometime discussed in terms 
of the standard thin Shakura-Sunyaev model. However, the standard thin disk 
is inconsistent with the super-Eddington accretion. Indeed, the radiative force 
cannot be greater than the effective gravity force, and at the surface of the 
standard disk this condition yields, 



F„ 



GMmp (H 
r 



>F r 



ad 



fx 



ad 



a T 3GMM 
c 8 7rr 3 



3r£ 
r 



(20) 



Indeed, all suggested "Kerr pseudo-potentials", are neither elegant, nor practical. 
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From i|2U|) Paczyhski ||j0] derived (H/r) max > M/Medd • This means that for a 
super-Eddington accretion with M/M^dd > 1 it also must be H/ r > 1. However, 
the standard model is a ID approximation to the 3D accretion physics and 
assumes H/r <§; 1. Its structure equations contain only zero and first order terms 
in H/r. Therefore, the super-Eddington accretion is outside of the standard 
model applicability range. Contrary to this, slim disks and Polish doughnuts are 
suitable to describe super-Eddington accretion: slim disks are accurate up to the 
second order terms in H/r, and Polish doughnuts are described by the full 3D 
equations (i.e. contain H/r terms of all orders). 

6.3 Adios: reasonable conclusions from false arguments 

The Blandford & Begelman "adios" paper [7^1 claims that accretion with small 
efficiency (Polish doughnuts, adafs) must necessarily experience strong outflows, 
because matter in these flows has everywhere a positive Bernoulli constant. 

While the very existence of strong outflows from small efficiency accretion 
black hole flows seems to be supported by observations, and therefore is most 
probably true, the reasons given by Blandford & Begelman to explain the out- 
flows, are certainly not correct, because, as discussed in [HU| by Abramowicz, 
Lasota & Igumenshchev, (1) A positive Bernoulli constant is only a necessary, 
but certainly not a sufficient, condition for outflows 16 . (2) The very argument 
that low efficiency accretion has a positive Bernoulli constant everywhere is not 
correct itself. It follows from an ad hoc mathematical assumption (made in order 
to make the problem easier to solve) that inefficient accretion flows are self- 
similar. Real flows obey inner and outer boundary conditions and for this reason 
are not self-similar. The boundary conditions imply negative Bernoulli constant 
at least close to both boundaries. 

7 Summary and conclusions 

Observations provide a clear evidence for super-Eddington luminosities powered 
by accretion onto black holes. Theory, solidly based on standard physics, predicts 
that such super-Eddington luminosities imply highly super-Eddington accretion 
rates and a very low efficiency of the black hole accretion. The powerful and 
simple Soltan argument shows, however, that the efficiency cannot be low all 
the time, because in the case of AGN this would be in a direct conflict with 
observations. 

Thus, observations and theory together seem to point that highly super- 
Eddington accretion, that really occur in several black hole sources, is a transient 
"flip-flop" phenomenon. A physical mechanism (or mechanisms ?) for the flip- 
flop behavior, in which periods of highly supper-Eddington but low efficiency 

16 The "a-positive-Bernoulli-implies-strong-outflows" is a wide spread fallacy. Rather 
surprisingly, as it is plainly incompatible with the well-known classic Bondi solu- 
tion, in which fluid has a positive Bernoulli constant everywhere, but experiences no 
outflows 69 . 
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accretion alternate with periods of highly efficient, sub-Eddington or ncarly- 
Eddington accretion is not yet known. Its understanding presents a great chal- 
lenge for all of us. 

Neither "adios", nor any other theoretical model, could at present explain 
strong outflows and jets that are observed in several supcr-Eddington black hole 
sources. Understanding of such outflows is another great challenge. 

In the context of the cosmic structure formation, one should conclude that 
there is no obvious reason for assuming that MEdd is the upper limit for the 
growth rate of the seed black holes. 
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